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Abstract. In this paper we present a method to solve 
the G-dwarf problem in the frame-work of analytical mod- 
els (based on the instantaneous recycling approximation, 
IRA). We consider a one-zone closed model without in- 
flows or outflows. We suppose a time-dependent Initial 
Mass Function (IMF) and we find an integral-differential 
equation which must be satisfied in order to honour the 
G-dwarf metallicity distribution as a function of the oxy- 
gen abundance. IMFs with one and two slopes are given 
and discussed also in the framework of a numerical chem- 
ical evolution model without IRA. We conclude that it is 
difficult to reproduce other observational constraints be- 
sides the G-dwarf distribution (such as [j^] vs [-7^]), and 
that an IMF with two slopes, with time-dependent shape 
at the low mass end, would be required. However, even in 
this case the predicted oxygen gradient along the disk is 
flat and radial flows would be required to reproduce the 
observed gradient. 
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1. Introduction 

It is well known that the metallicity distribution of G- 
dwarfs in the solar neighbourhood shows a deficit of 
metal-poor stars relative to the predictions of the Sim- 
ple Model of chemical evolution. This is the so-called G- 
dwarf problem, originally noted by van den Bergh (1962) 
and Schmidt (1963). 

Many solutions have been proposed, some in the frame- 
work of analytical models, based on the IRA. The problem 
is still present if we adopt the oxygen abundances in place 
of metallicity. 

An important result in the frame- work of IRA models 
with gas flows was pointed out by Edmunds (1990). He 
found that the G-dwarf problem cannot be solved by any 
outflow but it is possible to solve it with particular forms 
of inflow. Lyndon Bell (1975) 'best accretion model' and 
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Clayton's models (Clayton 1988) are models of just this 
kind. 

Of course it is also possible to reduce the number 
of metal-poor stars with respect to that predicted by 
the Simple Model by assuming metal dependent stellar 
yields (Maeder 1992) (more metals are produced at lower 
metallicity with a consequent Prompt Initial Enrichment 
(P.I.E.) (Truran and Cameron 1971)). However, it has 
been shown by several papers (Giovagnoli and Tosi, 1995; 
Carigi 1996), that in this way shallow gradients along the 
Galactic disk are produced, which do not agree with the 
most recent observational estimates (see Matteucci and 
Chiappini 1999 for a review). Actually, this is an unavoid- 
able problem still present in the models discussed in this 
paper where we consider time-dependent IMFs. 

It is still possible to obtain a Prompt Initial Enrich- 
ment by assuming an IMF variable with time and in par- 
ticular very flat at early times, in order to favour the for- 
mation of massive stars. 

In this paper we exploit the previous idea to investi- 
gate the time behaviour of the IMF. We assume the same 
hypothesis of the Simple Model (Tinsley 1980) with the 
exception of adopting a time-dependent IMF. In such a 
way it is possible to find an equation (Section 2.1) which 
must be satisfied in order to reproduce the observed distri- 
bution of G-dwarfs. We then use this equation to investi- 
gate the time behaviour of IMFs with one and two slopes. 
In Section 3 we use a numerical model without IRA to 
test the IMFs on other observational constraints. 

2. Recovering the history of the IMF 

In the following we need an analytical expression to ap- 
proximate the data concerning the metallicity distribution 
of G-dwarfs in the solar neighbourhood. The data are plot- 
ted in Fig. 1 where the differential distribution of oxygen 
abundances for the solar cylinder from Rocha-Pinto and 
Maciel (1996) is given (triangles). 

We have made use of the oxygen abundance rather 
than the iron abundance, in order to be consistent with 
chemical evolution models which use the instantaneous 
recycling approximation. Indeed, IRA is a good approxi- 
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Fig. 1. Metallicity distribution of 287 G-dwarf stars, from 
Rocha-Pinto and Maciel 1996 (triangles) and the function 
/($) (eq. 2) (continuous line). $ is the oxygen abundance 
in units of the solar abundance. AN* is the number of 
stars with oxygen abundance between <I> and $ + A$ 

mation for oxygen which is produced on short timescales 
(few million years) by supernovae (SNe) II as opposed to 
iron which is produced on long timescales (up to a Hubble 
time) by SNe la (Matteucci and Greggio 1986). In order 
to compare theory and data we adopt a simple relation 
between O and Fe reproducing the observed behaviour of 
these elements in the solar neighbourhood, the same as 
given by Pagel (1989): 
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The continuous line in Fig. 1 is the best function approx- 
imating data points according to the theory discussed in 
the appendix. 
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on the adopted method. Since the aim of this paper is to 
recover the history of the IMF from the data in Fig. 1, 
we need a method that does not introduce any a priori 
assumption on the specific functional form of the relation 
ln(AN*/ AQ). In other words our approach must be inde- 
pendent of any physical hypothesis on the star formation 
law at the basis of the observed data, and this is just the 
essence of the method discussed in the appendix. 

2.1. Basic Equations 

The basic assumptions of our models are the same of the 
Simple Model (Tinsley 1980) with the only exception of 
adopting a time-dependent IMF. Therefore we consider a 
model with the following properties: 

1. The system is one-zone and closed, namely there are 
no inflows or outflows 

2. The initial gas is primordial (no metals) 

3. The gas is well mixed at any time 

4. The instantaneous recycling approximation holds 
(namely the lifetimes of stars above 1M Q are negligible 
whereas those of stars below 1M Q are larger than the 
age of the Galaxy). 

5. The IMF is time-dependent, i.e. ip = <p(m, t), with the 
following normalization at any time 
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<p(m, t)mdm = 1 



(4) 



Under the previous assumptions the oxygen abundance $ 
(defined by eq.(l)) in the interstellar medium is governed 

by 



d<& = -^=-ads 

gZ 



(5) 



where Z = ~ 9.54 • 10 , g is the gas mass, 

pads = p'ds is the mass of oxygen produced and almost 
immediately returned into the interstellar medium when 
ds of interstellar material goes into stars. A fraction a 
of the mass gone into stars is not returned to the inter- 
stellar medium, but remains in long-lived stars or stellar 
remnants. It has also been assumed that 4>Z << 1. We 
have: 



rlviQ roa 

a = I ip(m,t)mdm+ / ip(m,t)m rem dm 
Jo Jm q 



Po 



a- = 0.454; $i = 0.690; $ 2 = 1-06; wl = 4.12; w2 = 2.61;(3) 

We have used the sophisticated method in the appendix to 
approximate the data because of the strong dependence of 
our results on the chosen function /($) and in particular 
on the value at $ = 0, /(0). Indeed, the point $ = is 
outside of the range of the available data and the metal- 
licity distribution extrapolated to $ = depends strongly 



and 
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(7) 



where p Q is the fraction (by mass) of newly produced and 
ejected oxygen by a star of mass m. We used for p a the 
expression given by Woosley and Weaver (1995) whereas 
for m rem the expression given by Tinsley (1980). 
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In our model the gas and stellar masses are related by 

dg = —ads (8) 

Now we want to relate the previous quantities to the ob- 
served metallicity distribution of G-dwarfs (/($) in eq. 
(2)) in order to find an equation for (p(m,t). We have: 



/(*) =ln(a 



where 



ds 



(9) 



a c = 



/■I.IMq 

= / Lp(m,t)mdm 

J0.8M Q 



(10) 



because the stars in our sample are in the range 0.8 — 
1.1M . 

Equations (5) and (9) give: 



/($) = lna c — Ina + Ing — Inp + InZo 



(11) 



On the other hand we can derive Ing from (5) and (8) in 
the following way: 



d<& 

Ing = Inga — Zq I — 
/o P 



(12) 



where go is the initial gas mass. Hence from (11) we obtain: 

/"* rf$ 

f ($) = lnot c — Ina — Z 1- lng a — Inp + lnZ (13) 

Jo P 

If we evaluate the previous equation at $ = (that is 
t = 0) we have: 



/(0) = lna c o + Ingo — lnp + InZo 



(14) 



where p' and a c o are the quantities in (6) and (10), re- 
spectively, with Lp(m,t) = ip(m,0). 

By differentiating eq. (13) with respect to $ we obtain 
finally 



— [lna c - Inp'} --Z = 

d$ P d$ 



(15) 



This is an integro-differential equation for the function 
<p(m, $) with the initial condition given by (14). The pre- 
vious problem has, of course, infinite solutions. Therefore 
to proceed further we have to make some assumption on 
the behaviour of the (p(m, $). In the next sections we shall 
investigate IMF with one slope (Sect. 2.2) and with two 
slopes (Sect. 2.3). 



2.2. Single Power-law IMF 

Let consider a single power-law IMF, namely 

^(m,$) = Cm-l 1+1 «l 



(16) 



where Ml and Mjj are respectively the smallest and the 
largest stellar mass (that we assume do not depend on the 
time) and the normalization is performed in the above 

ri l-x(<S>) 

mass range, i.e. C — 



Substituting eq. (16) in (15) we find the following equa- 
tion for x($): 



dx 



F 2 (x) + F 3 ($) 
F 1 (x) 



(17) 



which is a nonlinear I order differential equation with ini- 
tial condition given by substituting (16) in (14). The func- 
tions in (17) are: 



Fi (x) = -f- [lna c — Inp'] 
dx 
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Since F 2 (x) > 0, eq. (17) tells us that a constant IMF 
(i.e. with a slope x = x at any time) corresponds to a 
straight line for /($) with a negative slope (Indeed, in 
this case eq. (17) becomes F 3 ($) = % = ~F 2 (x )). This 
is the right behaviour since our model becomes the Simple 
Model when x = xo- 

The solutions x($) are plotted in Fig. 2 for two dif- 
ferent values of Mjj and Ml. We found very flat IMFs 
for low oxygen abundances (i.e. at initial times). At the 
solar metallicity ($ = 1) the slope is always steeper than 
a Salpctcr (1955) (x = 1.35). Decreasing Mjj decreases, 
of course, x, i.e. the IMF becomes flatter. The effect of a 
change in the Ml value is negligible (especially at low $). 

2.3. IMF with two slopes 

Let now consider the following IMF: 

TO -[i+*i] ifm<M 



<p(m, $) - C | MX2 _ Xlm _ [1+X2] iim>M 

The above IMF depends on the three parameters xi, x 2 
and M. We shall investigate the three following cases: 

1. M = M($), x 2 = 1.35 for values of xi ranging from 
— 1 to 0.2 (this IMF is similar to the one proposed by 
Larson (1998)); 

2. x\ — Xi($), x 2 = 1.35 for values of M ranging from 2 
to 10M Q ; 

3. x 2 = x 2 ($), for values of x\ ranging from —1 to 1. 
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Fig. 2. IMF slope (x($)) as a function of oxygen abun- 
dance (<&) for two choices of the parameters Mjj and Ml ■ 
The values of Mjj and Ml are given in solar mass. 
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Fig. 3. M($) in solar masses as a function of <f> for differ- 
ent x\. The slope x-i is always 1.35. 

In all the previous cases equation (15) gives us a non- 
linear I order differential equation for the ^-dependent pa- 
rameter, with the initial conditions given by the equation 
(14). We consider always M v = 100M Q and M L = 0.1M© 
since the dependence on these parameters is negligible. 
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Fig. 4. x\ ($) as a function of $ for different M. x 2 — 1.35. 

The function M(<&) in the first case considered is 
shown in Fig. 3. At the initial time (<j> = 0) the mass 
M increases by increasing the slope at low mass end. In 
particular when x\ = 0.2 M(0) ~ 100M Q = M v and 
therefore there is no solution for x\ greater than 0.2. 

Fig. 4 shows the slope at the low mass end (xi($)) 
for values of M in the range 2 — 10M Q . The initial values 
of x\ are always negative and moreover they increase by 
increasing M, as expected. 

Finally Fig. 5 shows the results for the third case con- 
sidered. Here the mass M is chosen in order to reproduce 
a final slope x 2 — 1-35. 

3. Application of the derived IMFs to a numerical 
model 

Both the proposed IMFs can reproduce the observed G- 
dwarf distribution and therefore we want to test the va- 
lidity of the inferred IMF by using numerical models of 
galactic chemical evolution . This can be done by study- 
ing the effect of the above IMFs on the chemical evolution 
of the solar neighbourhood. The model used is that of 
Matteucci and Francois (1989) where a detailed descrip- 
tion can be found. The main difference with that model is 
that we assume here a very rapid formation of either the 
halo and the disk thus simulating a closed model. This is 
required by the fact that we derived the IMF under the 
assumption of a closed model. In the original model of 
Matteucci and Francois (1989) the timescale for the for- 
mation of the solar neighbourhood was about 3 — 4 Gyr 
and it was chosen in order to best fit the G-dwarf metal- 
licity distribution of Pagel and Patchett (1975), under the 
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Fig. 5. iE2( < E > ) as a function of <f> for different x\. M is 
chosen in order to have a final slope X2 — 1-35 

assumption of a constant IMF. Recently, Chiappini et al. 
(1997) presented a more realistic model for the Galaxy 
where the evolution of the halo-thick disk and the thin 
disk are decoupled in the sense that the halo-thick disk is 
formed on a short time scale of the order of 1 — 2 Gyr, 
whereas the thin disk is assumed to form very slowly. 

In particular, in order to best fit the G-dwarf distri- 
bution of Rocha-Pinto and Maciel (1996), Chiappini et al. 
(1997) found that, with an IMF constant in time, a time 
scale of about 8 Gyrs is required for the formation of the 
solar neighbourhood. It is interesting to check if the de- 
rived IMFs are able to reproduce other observational con- 
straints besides the G-dwarf distribution, for example the 
[O/Fe] vs [Fe/H] trend which is normally very well repro- 
duced by models with constant IMF and infall and taking 
into account detailed nucleosynthesis from type II and la 
SNe. In this framework, in fact, the plateau shown by the 
data for [Fe/H] < —1.0 is interpreted as due to the pol- 
lution by type II SNe which produce an almost constant 
[O/Fe] ratio. The subsequent decrease of the [O/Fe] ratio 
for [Fe/H] > —1.0 is then due to the occurrence of type 
la SNe exploding with a temporal delay relative to SNe 
II. In the Matteucci and Francois (1989) and Chiappini 
et al. (1997) model the SNe la are supposed to originate 
from white dwarfs in binary systems following the formal- 
ism of Matteucci and Greggio (1986). In Fig. 6 we show 
the G-dwarf metallicity distribution as predicted by the 
numerical model and, as expected, the agreement is quite 
good. However, in Fig. 7 we show the predicted [O/Fe] 
versus [Fe/H] relation and the agreement with the ob- 
servations is poor, especially in the domain of metal poor 




| Fe/H] 

Fig. 6. G-dwarf metallicity distribution as predicted by 
the numerical model (solid line) with the single power law 
IMF discussed in Sect. 2.1, and from experimental data 
(Rocha-Pinto and Maciel 1996) 



stars where the new IMF predicts an increasing amount 
of massive stars. On the other hand, the slope for disk 
stars is well reproduced. Another problem with this IMF 
is the high O and Fe content reached by the model at 
the sun birth and at the present time. These high abun- 
dances are not evident in Fig. 7 since the abundances are 
normalized to the predicted solar abundances. Finally, an- 
other problem is that Fig. 7 shows that [Fe/H] starts de- 
creasing more than oxygen after having reached the solar 
abundance. This is due to the large dilution from dying 
low mass stars and to the fact that the Fe abundance de- 
creases more rapidly than that of oxygen, due to the fact 
that O is continuously produced, although at a low level, 
by SNe II on very short timescales. Iron, on the other 
hand, comes mostly from type la SNe born at early times 
when the IMF was top-heavy favoring massive stars rel- 
ative to the type la SN progenitors. The predicted solar 
abundances, namely the abundances at 4.5 Gyrs ago, are 
X = 2.6 • 10~ 2 and X Fe = 6.3 • 10~ 3 to be compared with 
the same abundances from Anders and Grevesse (1989) : 
X = 9.59 • 10~ 3 and X Fe = 1.17 • 10~ 3 (mass fractions). 

Concerning the IMF with two slopes (discussed in Sect. 
2.3), the G-dwarf metallicity distribution as predicted by 
the numerical model is again in good agreement with the 
observations, in all the cases considered. On the other 
hand, the predicted [O/Fe] versus [Fe/H] is still in poor 
agreement with the observations with the exception of the 
second case (x\ — #i(3>)) for which we give the G-dwarf 
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Fig. 7. [O/Fe] vs [Fe/H] as predicted by the numerical 
model (solid line) with the single power law IMF (Sect. 
2.1), and as observed 

metallicity distribution and the [O/Fe] versus [Fe/H] in 
Figs 8 and 9, respectively, for M — 5M Q (we found results 
very similar to the ones shown in the Figs 8 and 9 for val- 
ues of M in the range 2 — 10M©). This time-dependent 
IMF fits the [O/Fe] vs [Fe/H] relation much better than 
the other time dependent IMFs, since the variation in the 
number of type II SNe at early times is less than in the 
other cases. However, it is clear that the best agreement 
with the data (especially for [Fe/H] < —1.0) is achieved 
with a constant IMF. 

We also investigated the case of an IMF with only the 
lower mass limit dependent on <j> ((Ml = Ml($)). How- 
ever, in this case, the equations (14) and (15) give solu- 
tions which are not physical. In fact, in order to lower 
the number of stars in the range of mass 0.8 — 1.1M Q at 
early times, to solve the G-dwarf problem, our equations 
fix the value of M L {$) (for $ ~ 0) very close to 1.1M©. 
This result does not have any physical meaning because it 
is strongly dependent on the particular mass range of in- 
terest. This is the reason why the inferred solution, when 
used in a numerical model of chemical evolution, it is not 
able to reproduce the observational G-dwarf distribution. 

The predicted solar abundances related to the IMF 
adopted in Figs 8 and 9 are still larger than the observed 
ones (X = 3.1 • 10~ 2 and X Fe = 8.6 • 10~ 3 ). However, 
these absolute values are strongly model dependent and 
in particular it is possible to lower these abundances by 
increasing the infall. 



Fig. 8. G-dwarf metallicity distribution as predicted by 
the numerical model (solid line) and from experimental 
data (Rocha-Pinto and Maciel 1996). The IMF adopted is 
with two slopes (discussed in Sect. 2. 2) with x\ = xi(<j>) 
x 2 = 1.35 and M = 5M Q 

4. Conclusions 

We proposed a method to solve the G-dwarf problem in 
a closed box model with a time-dependent IMF, based on 
the IRA. The method gives us an equation (eq. 15) which 
has infinite solutions. Therefore, in order to use this equa- 
tion, we had to make some assumptions on the behaviour 
of the IMF. In particular, we considered both a single 
power-law IMF and an IMF with two slopes. 

We tested the validity of the inferred IMF by using 
numerical models of galactic chemical evolution, namely 
by relaxing the IRA, and we find the following results: 

1. all the IMFs investigated can reproduce the observed 
G-dwarf distribution; 

2. a single power-law IMF fails in reproducing the be- 
haviour of abundances (in particular the [O/Fe] vs 
[Fe/H] relation); 

3. in order to reproduce the behaviour of abundances be- 
sides the G-dwarf problem, an IMF with a time depen- 
dence at the low mass end is required. However, the fit 
produced by such an IMF of the [O/Fe] vs [Fe/H], is 
not as good as that produced by a constant IMF. More- 
over, this IMF, like most of the variable IMF proposed 
insofar, fails in reproducing the oxygen gradient along 
the Galactic disk, unless other assumptions such as in- 
creasing star formation efficiency with galactocentric 
distance and/or radial flows are introduced. We com- 
puted the expected gradients of oxygen along the disk 
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Fig. 9. [O/Fe] vs [Fe/H] as predicted by the numerical 
model (solid line) and as observed. The adopted IMF is 
the same as in Fig. 8. 

by adopting this IMF and the model of Chiappini et 
al. (1997). We found that the O gradient disappears. 
A variable efficiency of star formation could recover 
the gradient but at expenses of the gas distribution 
along the disk. We did not try to include radial flows 
but the conclusion that a variable IMF of this kind 
worsen the agreement with the disk propeties relative 
to a constant IMF seems unavoidable. 

4. an IMF similar to the one proposed recently by Lar- 
son (1998), fails in reproducing the [O/Fe] vs [Fe/H] 
relation because of the strong time dependence of the 
number of type II SNe at early times. 

5. In summary, from the analysis of models with variable 
IMF done in this paper, we are tempted to conclude 
that infall models with a constant IMF are the best 
solution of the G-dwarf problem, since variable IMFs 
could in principle solve it but they are not able to 
reproduce the properties of the galactic disk. 
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5. Appendix 

We give a method to approximate observational data in 
Fig. 1, based on the regularization theory of Tikonhov 
(1963). 

The method is general and it applies whenever one 
wants to approximate some data by an analytical function, 
and nothing is known about the physics at the basis of 
these data. 

Accordingly to Tikonhov's regularization theory 
(1963) the function /(<&) which approximates the data 
in Fig. 1 is determined by minimizing a cost functional 
E[f], so-called because it maps functions (in some suit- 
able function space) to the real line. E[f] is the sum of 
two terms 



E[f]=E s [f] + E c [f] 



(21) 



where E s [f] is the standard error term and E c [f] the reg- 
ularizing term. The first-one measures the standard error 
(distance) between the desired response fi and the actual 
response f($i), 



1 N 

E s [.f] = 7 :Y.^-f^)) 2 



(22) 
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where N is the total number of available data (in our 
case N = 12). The second-one depends on the geometric 
properties of the approximating function /(<&), 

E c [f) = \\\Pff (23) 

where P is a differential operator. As suggested by Pog- 
gio and Girosi (1990) the best choice for P consists in a 
differential operator invariant under both rotations and 
translations. This is defined by 

||P/|| 2 = f>||£>*/|| 2 (24) 

fc=0 

where cik = |^ with a a constant associated with the 
data point $^ and 

l|DW = /I(fp)* (25) 

The function which solves the variational problem given 
by eq. (21) (when the regularizing term is specified by eq. 
(24)) is: 

/(*) = 5^- (26) 

3=1 

which consists of a linear superposition of multivariate 
Gaussian basis functions with centers $j . The above the- 
ory establishes only the form of the function /(<£>), but 
does not solve completely our problem. We do not know 
for example the number M of gaussians in (25). 

As suggested by Guyon et al. (1992) and Vapnik 
(1992), the method to fix the parameters in the expres- 
sion (26) is to minimize the sum of a pair of competing 
terms. The former is again the standard error given by eq. 
(22), and decreases monotonically as the number of pa- 
rameters is increased. The latter measures the complexity 
of the model and increases by increasing the number of pa- 
rameters. Therefore there is an optimal compromise which 
minimizes the sum. It is possible to demonstrate (Amari 
et al. 1997) that we obtain this compromise by using N — k 
data and the average error made on the patterns left out. 
This method is called leave-k-out cross validation. When 
N is small, the most reasonable choice is k = 1 (leave-one- 
out). 

Since in our case the value of N is small we apply the 
leave-one-out cross-validation method in order to fix the 
parameters in (26). The average errors made on the left 
data point are reported in Table 1 related to several values 
of M. 

Consequently the best model occurs for M = 2 with 
values of the parameters appearing in eq. (26) given by 
eq. (3). 



M 


Error 


1 


0.9131 


2 


0.4988 


3 


0.7545 


4 


0.8080 


5 


0.9097 



Table 1. Average errors made on the left data points fol- 
lowing the leave-one-out cross validation method for sev- 
eral values of M. 



